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INTRODUCTION 

Nondestructive  evaluation  (NDE)  of  small  spaces  such  as  internal  holes  or  pipes  in 
engineering  structures  is  generally  a  challenging  task.  The  difficulty  usually  arises  due  to  the  small 
space  to  insert  a  probe  and  the  consequent  impediments  to  manipulate  the  probe  in  the  tight  space 
and  acquire  reliable  data  at  high  speed  with  precision  and  accuracy.  Laser-based  inspection 
systems  have  been  reported  for  inspection  of  pipes  of  both  large  [~20  to  25  cm  in  diameter  (refs.  1 
and  2)]  and  small  [<1  cm  in  diameter  (ref.  3)].  In  addition,  a  computer  vision-based  inspection  system 
(ref.  4)  and  a  three-dimensional  vision  inspection  using  structured  light  (ref.  5)  have  also  been 
proposed  for  inspecting  internal  surfaces  of  holes  that  have  a  diameter  of  <1  cm.  Notwithstanding 
these  advances,  the  inspection  of  small  spaces  is  further  exacerbated  if  other  surface  imperfections 
such  as  stress  corrosion  pits  are  also  present.  In  these  circumstances,  detection  of  surface  flaws, 
such  as  cracks  in  small  internal  surfaces  often  initiated  at  corrosion  pits  (ref.  6),  is  almost  impossible. 
Thus,  unambiguous  identification  of  surface  flaws  in  the  midst  of  corrosion  pits  and  the  derivation  of 
accurate  and  reliable  statistical  parameters  of  their  distribution  are  an  important  problem  for  routine 
NDE  of  small  spaces.  Herein,  the  application  of  a  modified  commercial-off-the-shelf  (COTS)  laser- 
based  inspection  system  (ref.  3)  to  probe  the  internal  surface  textures  of  a  <1-cm  diameter  hole 
inside  a  steel  structure  is  described.  Specifically,  signal  processing  and  data  analysis  techniques 
have  been  developed  so  as  to  derive  reliable  statistical  parameters  of  the  inside  surface  textures  of 
the  holes  after  correcting  for  laser  probe  wander  and  deviations  in  the  geometry  of  the  hole. 


EXPERIMENTAL  PROCEDURE 

The  holes  that  were  examined  were  of  ~6  mm  in  diameter  and  ~1 9  mm  in  length  embedded 
in  a  tubular  section  made  of  high  strength  steel.  The  surface  texture  of  the  inside  of  the  holes  was 
probed  using  a  modified  COTS  laser-based  probe,  the  BEMIS  SC™  Small  Caliber  Inspection  System 
of  Laser  Techniques  Company,  Redmond,  WA  (ref.  3),  as  shown  in  figure  1 .  The  tubular  sections 
were  positioned  on  X-Y  stages  for  easy  manipulation.  The  laser  probe  was  aligned  normal  to  the 
hole  and  centered  as  much  as  possible,  and  the  start  of  the  scans  was  at  a  12  o'clock  position  of 
each  hole.  The  linear  and  rotary  resolutions  were  set  at  50.8  and  25.4  pm,  respectively.  By 
experimentation,  it  was  found  that  reproducible  results  were  obtained  when  the  laser  power  was  set 
in  auto  mode  and  the  detector  gain  was  set  at  5.  The  raw  data  files  were  imported  to  MATLAB  for 
further  signal  processing.  Some  of  the  tubular  sections  were  subjected  to  a  safety  and  reliability  test 
simulating  storage  conditions  to  withstand  exposure  to  repeated  cycles  of  extreme  temperatures  and 
humidity  as  per  MIL-STD-331C  test  Cl  (ref.  7).  The  tubular  sections  were  exposed  to  a  28-day  (two 
14-day  cycles)  schedule  of  temperature  and  humidity  variations  to  induce  corrosion  in  the  steel.  The 
temperature  extremes  were  +71°  and  -54°C,  and  a  relative  humidity  of  95%  was  used  at  +71  °C.  As 
a  rule  of  thumb,  this  testing  in  an  unpackaged  configuration  is  considered  equivalent  to  2  years  of 
storage  in  an  uncontrolled  environment  (ref.  8).  The  tubular  sections  were  subjected  to  a  total  of 
three  28-day  cycles  for  an  equivalent  exposure  of  6  years  in  an  uncontrolled  environment  and  were 
scanned  after  28,  56,  and  84  days.  Following  the  treatment,  the  corrosion  build-up  inside  the  holes 
was  removed  using  brass  bore  brushes  and  ultrasonic  cleaning.  Due  to  the  increased  corrosion 
residue  buildup  after  the  84-day  treatment,  stainless  steel  brushes  in  place  of  brass  and  a  more 
rigorous  cleaning  cycle  were  employed  compared  to  the  28  and  56-day  treatments.  This  cleaning 
method  was  repeated  as  needed  such  that  most  of  the  corrosion  residue  was  removed  and  the  scan 
results  did  not  vary  appreciably  after  successive  cleaning  steps. 
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Figure  1 

Experimental  setup  of  probing  the  inside  surfaces  of  holes  in  a  steel  tubular  section  using  a  modified 

COTS  laser-based  probe 


RESULTS  AND  DISCUSSION 

Prior  to  extracting  statistical  information  of  surface  features,  the  following  operations  were 
first  applied  to  the  raw  scan  data: 

•  Each  scan  consisted  of  a  raster  of  measured  radii  r(z,  0)  as  a  function  of  linear  dimension 
(z)  and  angle  (0).  Minimum  and  maximum  values  for  linear  dimension  (z)  were  selected, 
and  the  scan  is  clipped  to  these  bounds  to  limit  the  analysis  to  points  inside  the  hole.  While 
there  is  some  inherent  subjectivity  in  this  task,  every  attempt  was  made  to  remain  as 
consistent  as  possible.  The  plot  in  figure  2  illustrates  the  chosen  limits  for  a  particular  scan. 
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Figure  2 

Plot  illustrating  the  chosen  limits  (the  dotted  lines)  for  a  particular  scan 

After  clipping,  each  scan  was  median  filtered  by  applying  a  two-dimensional  median  filter. 
That  is,  each  pixel  r[n,  m]  is  replaced  with  the  median  value  of  the  following  set  of  points 
{r[n,  m],  r[n-1,  m],  r[n+1,  m],  r[n,  m-1],  r[n,  m+1]}.  This  operation  is  quite  effective  in 
removing  the  speckle  noise  (i.e.,  isolated  single  point  outliers)  that  is  most  likely  associated 
with  the  measurement  instrument. 

A  correction  was  applied  to  the  measured  radii  to  correct  for  the  probe  deviation  off  the 
center  axis  of  the  hole.  At  each  slice  in  z,  an  ellipse  was  fit  to  the  measured  radii.  Each 
radius  measurement  was  then  replaced  with  the  distance  to  the  center  of  the  ellipse,  which 
was  taken  as  the  true  axis  of  rotation.  The  parameters  of  this  ellipse  were  estimated  by 
minimizing  the  sum  of  the  squared  algebraic  distances.  If  the  parameters  of  the  ellipse  are 
denoted  as  a  =  [a,  b,c,d,e,f]  and  the  sampled  points  by  {p1(p2,  --- Pjv)  with  p„  =  (*n<yn)> 
then  the  parameters  are  estimated  by  finding  the  coefficients  that  minimize  the  following: 


I£=i  F{a,pn? 


(1) 


where  F(a,pn)  =  ax %  +  bxnyn  +  cy%  +  dxn  +  eyn  +  f  =  0.  This  method  appears  to  be 
effective  in  correcting  the  meandering  of  the  probe  as  it  goes  through  the  hole.  For 
example,  shown  in  figure  3a,  is  a  scan  that  illustrates  the  effects  of  a  slight  probe  wander. 
Because  the  probe  is  slightly  off  the  center  axis  of  the  hole,  the  measured  radius  is  larger  at 
certain  angles  and  smaller  at  others.  After  applying  the  ellipse  fitting  technique,  the 
resulting  scan  is  much  more  cylindrical  while  still  retaining  much  of  the  small-scale  features 
as  shown  in  figure  3b. 
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(a) 

Scan  showing  effects  of 
probe  wander  as  the  probe 
is  slightly  off  the  center 
axis  of  the  hole 


(b) 

Scan  illustrating  the 
efficacy  of  the  ellipse  fitting 
technique  in  removing 
the  probe  wander 


Figure  3 

Effects  of  probe  wander  and  ellipse  fitting  technique  in  removing  probe  wander 

After  correcting  for  the  probe  wander,  the  angular  values  (0)  were  replaced  with  an 
“unwrapped”  linear  dimension  (x)  using  a  nominal  radius  using  the  mapping  x  =  rref9, 
[where  vref  was  taken  as  2.78  mm  (nominal  reference  radius),  and  9  is  the  measurement 
angle  in  radians]. 

The  scan  was  resampled  (using  bicubic  interpolation)  for  convenience  to  a  common 
resolution  in  both  x  and  z. 

It  was  also  useful  to  filter  out  lower  frequency  components.  A  two-dimensional,  circularly 
symmetric,  low-pass  filter  was  applied  to  the  measurements  to  generate  a  background 
image  consisting  of  the  slowly  varying  components.  The  zero  frequency  component  was 
retained  in  order  to  put  the  resulting  signal  on  the  same  scale  as  the  original  (i.e.,  around 
2.78  mm).  This  was  accomplished  by  removing  all  spatial  frequencies  less  than  0.25  mm~ 1 
in  all  directions  via  spatial  filtering.  Since  potential  pits  are  expected  to  be  highly  localized  in 
space  (high  spatial  frequency)  and  because  of  the  possibility  of  longer  wavelength 
variations  (e.g.,  deviations  in  the  geometry  of  the  hole,  remaining  probe  wander),  this 
background  provided  a  convenient  baseline,  removing  any  slowly  varying  spatial 
components  to  compare  the  measured  radii  for  detection  of  pits. 

A  binary  image  was  formed  from  the  measured  radii  that  are  at  a  given  threshold  above  the 
background  image. 
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•  A  single  pixel  above  the  background  is  most  likely  not  a  pit.  To  remove  these  potential  false 
positives,  a  morphological  “opening”  operation  is  performed  to  the  previously  generated 
binary  image.  This  consisted  of  “erosion”  (retaining  only  points  such  that  all  points  in  a 
stencil  are  present)  followed  by  “dilation”  (retaining  points  with  any  point  in  a  stencil  being 
present).  What  this  operation  accomplishes  is  to  remove  groups  of  points  below  a  given 
radius. 

•  The  remaining  points  in  the  binary  image  are  then  associated  to  common  pits.  In  order  to 
associate  pixels  together,  a  criteria  for  connectedness  is  required.  The  criteria  chosen  was 
8-connected  (i.e.,  a  pixel  r[n,  m]  can  be  connected  to  a  set  of  pixels  if  any  of  its  eight 
surrounding  neighbors  are  a  member  of  the  set). 

Figure  4  illustrates  successively  processed  plots  of  the  original  measurement  after  steps  6 
(fig.  4a),  7  (fig.  4b),  8  (fig.  4c),  and  9  (fig.  4d).  Figure  4a  shows  the  successively  processed  plots  after 
background  image  determination  by  removing  all  spatial  frequencies  less  than  0.25  mm-1  in  all 
directions  via  fast  Fourier  transform  filtering.  Figure  4b  shows  the  successively  processed  plots 
forming  a  binary  image  from  the  measured  radii  that  are  above  a  threshold  of  0.125  mm.  Figure  4c 
shows  the  successively  processed  plots  after  a  morphological  operation  of  the  binary  image  shown 
in  figure  4b.  Figure  4d  shows  the  association  of  points  to  common  pits  using  the  8-connected  criteria. 


Background  (Cutoff:  0.251mm'1]) 
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Figure  4 

Successively  processed  plots  of  the  original  measurement 


Scan  Statistics 

The  effects  of  the  corrosion  treatments  on  the  inside  surfaces  of  the  holes  can  be  best 
captured  using  the  statistical  information  of  various  surface  texture  parameters  as  defined  in  the 
International  Organization  for  Standardization  (ISO)  251 78  for  quantification  of  areal  surface  texture 
(refs.  9  and  1 0).  The  cumulative  distribution  of  surface  heights  as  a  function  of  percentile  can  be 
illustrated  via  the  Abbott-Firestone  curves  (refs.  9  and  10)  after  applying  all  of  the  previously 
mentioned  operations  of  clipping,  median  filtering,  correcting  for  probe  wander,  unwrapping,  and 
removing  low  frequency  components  have  been  applied.  Figure  5  shows  the  Abbott-Firestone 
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curves  for  a  hole  after  corrosion  treatment  for  28,  56,  and  84  days  in  comparison  to  the  initial 
condition  illustrating  the  increase  in  surface  heights  with  the  number  of  days  of  treatment.  It  can  be 
seen  that  at  lower  percentiles  (e.g.,  5th  through  10th),  the  surface  height  tends  to  increase  with  time 
illustrating  the  likelihood  of  the  formation  of  pits.  It  was  observed  that  the  maximum  surface  height, 
however,  is  not  necessarily  always  after  the  longest  (i.e.,  84  days)  treatment  suggesting  perhaps  that 
there  is  material  buildup  covering  large  peaks  and  that  the  post-corrosion  cleaning  process  is  not 
adequate  in  removing  the  corrosion  residue  completely.  However,  there  is  a  clear  widening  of  the 
distribution  of  material  with  respect  to  the  number  of  days  of  treatment,  exhibiting  both  large  heights 
(potential  pits)  as  well  as  low  heights  (buildup  of  material)  at  longer  times. 


Abbott-Firestone  Curve 


Percentile 


Figure  5 

Abbott-Firestone  curves  of  a  hole  that  is  subjected  to  corrosion  treatment  for  28,  56,  and  84  days 

Pit  Definition 

Identification  of  pits  inside  the  holes,  especially  after  the  corrosion  treatment,  is  an  important 
engineering  problem  in  assessing  the  engineering  life  of  structures  as  most  often  surface  flaws  such 
as  cracks  are  initiated  at  corrosion  pits  (ref.  6).  There  is,  however,  a  great  deal  of  difficulty  in 
ascertaining  what  constitutes  a  “pit”  in  a  given  laser  scan  of  the  inside  surface  features  of  a  hole, 
particularly  when  the  pit  is  located  deep  inside  the  hole  and  not  readily  available  for  visible 
inspection.  Clearly,  a  pit  should  be  localized  in  space  and  should  have  a  sufficient  height.  The 
requirement  that  pits  be  localized  in  space  is  one  of  the  primary  reasons  for  removing  lower 
frequencies  (i.e.,  spatial  frequencies  less  than  0.25  mm'1). 

There  are  several  parameters  that  affect  the  number  of  pits  detected.  The  most  critical  of 
which  are:  (1)  detection  threshold  value,  (2)  size  of  morphological  “opening”  stencil,  and  (3)  criteria 
for  connectedness.  In  order  for  a  pixel  to  be  classified  as  a  potential  pit,  the  pixel’s  value  must 
exceed  the  background  image  by  a  value  greater  than  or  equal  to  the  detection  threshold.  To  get 
realistic  values  for  the  detection  threshold,  a  tubular  section  that  had  been  exposed  to  the 
environment  for  several  years  was  cut  open  and  the  pit  depths  were  obtained  using  an  optical 
comparator.  An  inspection  putty  was  used  to  form  a  mold  of  the  cut  opened  hole  surface,  and  after 
~15  min,  the  surface  features  were  impressed  on  the  putty.  The  optical  comparator  scan  of  the  putty 
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surface  indicated  that  the  deepest  pit  is  of  the  order  of  0.  675  mm  and  the  smallest  pit  is  of  the  order 
of  0.125  mm.  The  laser  scan  results  from  this  section,  obtained  by  using  a  threshold  of  0.125  mm, 
spatial  filter  of  0.25  mm-1,  and  a  morphological  opening  stencil  of  1  pixel  radius,  showed  that  the 
maximum  and  minimum  pit  depths  were  0.683  and  0.131  mm,  respectively.  The  maximum  value  is 
remarkably  close  to  that  obtained  by  the  optical  comparator  demonstrating  high  accuracy  of  the  laser 
scan  method  in  pit  depth  determination.  The  pit  distributions  were  computed  for  all  the  tubular 
sections  tested  in  the  study  setting  thresholds  at  -0.075,  0.125,  and  -0.25  mm.  For  the  results 
presented  in  figure  6,  unless  otherwise  stated,  the  following  parameters  were  held  constant: 


Bandwidth  of  low-pass  filter:  0.25  mm-1 
Morphological  opening  stencil  used  (radius  1  pixel): 


010 

111 

LoioJ 


Connectedness:  8-connected 


Pit  Depth  [um]  Pit  Depth  [um] 


Pit  Depth  [um] 
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Figure  6 

Binary  images  showing  the  progression  of  pit  formation  with  the  number  of  days  of  treatment  for  the 

same  hole  in  a  tubular  section 

The  progression  of  the  pit  formation  as  a  function  of  the  number  of  days  of  treatment  is 
illustrated  in  figure  6.  The  increased  density  of  pits  greater  than  0.125  mm  (threshold)  after  the  84- 
day  treatment  is  clearly  seen  in  the  figure.  The  top  left  image  is  the  initial,  the  top  right  is  after  28 
days,  the  bottom  left  is  after  56  days,  and  the  bottom  right  is  after  84  days.  The  increased  pit  density 
after  84  days  treatment  is  clearly  noticed  in  the  bottom  right  image.  The  numbers  on  the  images 
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denote  the  pit  depths.  For  these  images,  the  threshold  was  set  at  0.125  mm  and  the  morphological 
opening  stencil  of  1-pixel  radius  was  used. 

The  effects  of  the  threshold  and  the  morphological  opening  stencil  radius  on  the  pit 
distribution  were  examined.  Figure  7  illustrates  the  effect  of  the  threshold  value  on  the  ensuing  pit 
distribution  inside  a  hole  of  a  tubular  section  after  the  84-day  treatment  applied  to  the  same  scan 
data.  Points  that  exceed  the  background  protrude  through  the  surface  as  illustrated  in  figure  7,  and 
the  effect  of  the  three  different  thresholds  (0.25,  0.1 25,  and  0.075  mm)  on  the  ensuing  pit  distribution 
is  clearly  seen.  Figure  8  illustrates  the  effect  of  the  varying  threshold  for  a  fixed  morphological 
opening  stencil  radius  (1  pixel)  and  varying  morphological  opening  stencil  radius  for  a  fixed  threshold 
(0.1 25  mm)  on  the  number  of  pits  for  the  scan  data  from  one  of  the  tubular  sections  after  the  84-day 
treatment.  The  top  plot  in  figure  8  illustrates  the  effect  of  threshold  on  the  number  of  pits  for  a 
morphological  opening  stencil  of  1 -pixel  radius.  The  bottom  plot  shows  the  effect  of  varying  the 
morphological  opening  radius  on  the  number  of  pits  for  a  threshold  value  of  0.125  mm.  Both  the  plots 
were  derived  from  the  same  scan  data  from  one  of  the  tubular  sections  in  the  original  condition  and 
after  the  28,  56,  and  84-day  treatments.  The  progressive  increase  in  the  number  of  pits  with 
increasing  number  of  treatment  days  clearly  evince  the  increased  corrosion. 
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Figure  7 

Pit  distribution  inside  a  hole  of  a  tubular  section  after  the  84-day  treatment  for  three  different 

thresholds 


Approved  for  public  release;  distribution  is  unlimited. 

UNCLASSIFIED 

8 


UNCLASSIFIED 


(each  pixel  is  0.025mm  x  0.025mm;  Threshold  is  0.125 


Figure  8 

Effect  of  the  varying  threshold  for  a  fixed  morphological  opening  stencil  radius  (1  pixel)  and  varying 
morphological  opening  stencil  radius  for  a  fixed  threshold  (0.125  mm) 

Finally,  figure  9  illustrates  the  pit  density  as  a  function  of  pit  depth  for  the  different  number  of 
treatment  days  derived  from  the  data  obtained  from  all  15  of  the  tubular  sections  examined  in  the 
study.  The  number  of  pits  for  a  morphological  opening  stencil  of  1 -pixel  radius  is  calculated  in  25-pm 
bins,  and  this  is  divided  by  the  total  area  to  compute  the  pit  density  in  each  bin.  The  number  of  pits  in 
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the  25-pm  bins  is  shown  at  the  bottom  of  the  figure.  A  few  observations  can  be  made:  (1 )  pit  density 
decreases  monotonically  with  increasing  pit  depth,  (2)  the  maximum  pit  depth  beyond  which  the  pit 
density  is  negligible  increases  with  the  number  of  treatment  days,  and  (3)  for  a  given  pit  depth  below 
the  maximum,  the  pit  density  increases  with  the  number  of  treatment  days.  For  example,  for  a  pit 
depth  of  0.1  mm,  the  pit  density  increases  by  ~2.5  and  ~4.5  times,  respectively,  after  56  and  84  days 
compared  to  that  after  28  days  of  treatment.  Given  that  the  environment  test  of  temperature  and 
humidity  exposure  used  in  the  study  is  equivalent  to  ~2  years  of  uncontrolled  environment,  it  may  be 
surmised  that  for  every  two  years  of  exposure,  the  growth  rate  of  pits  is  approximately  doubled. 


100  150  200  250  300  350 


Pit  Depth  [pm] 


Figure  9 

Plot  of  pit  density  as  a  function  of  pit  depth  for  different  number  of  treatment  days 


CONCLUSIONS 

A  modified  commercial  off-the-shelf  laser  based  probe  was  applied,  the  BEMIS  SC™  Small 
Caliber  Inspection  System  of  Laser  Techniques  Company,  Redmond,  WA,  to  nondestructively 
examine  the  inside  surface  texture  of  holes  ~6  mm  in  diameter  and  ~19  mm  in  length  embedded  in  a 
tubular  section  made  of  high  strength  steel.  The  data  was  analyzed  using  a  signal  processing 
algorithm  developed  in  MATLAB  that  included:  (1)  correction  for  the  laser  probe  wander  of  the  center 
axis  of  the  hole,  (2)  a  low-pass  filter  to  remove  any  slowly  varying  spatial  components  and  provide  a 
convenient  reference  surface,  (3)  determination  of  surface  pits  that  are  above  a  give  threshold  value, 
(4)  a  morphological  operation  to  remove  potential  false  positives,  and  (5)  a  criteria  of  connectedness 
to  associate  pixels  to  common  pits.  The  tubular  sections  were  subjected  to  a  28-day  (two  14-day 
cycles)  schedule  of  temperature  and  humidity  variations  to  induce  corrosion  in  the  steel.  The 
temperature  extremes  were  +71°  and  -54°C,  and  a  relative  humidity  of  95%  was  used  at  +71  °C. 

The  cumulative  distribution  of  surface  heights  as  a  function  of  percentile,  illustrated  via  the  Abbott- 
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Firestone  curves,  revealed  that  the  surface  height  tends  to  increase  with  the  number  of  days  of 
treatment  suggesting  the  formation  of  pits.  Further  analysis  of  the  pit  distributions  showed  that  the  pit 
density  decreases  monotonically  with  increasing  pit  depth,  the  maximum  pit  depth  beyond  which  the 
pit  density  is  negligible  increases  with  the  number  of  treatment  days,  and  for  a  given  pit  depth  below 
the  maximum,  the  pit  density  increases  with  the  number  of  treatment  days.  For  example,  for  a  pit 
depth  of  0.1  mm,  the  pit  density  increases  by  ~2.5  and  ~4.5  times,  respectively,  after  56  and  84  days 
compared  to  after  28  days  of  treatment.  No  surface  defects  such  as  cracks  were  observed  in  any  of 
the  sections  even  after  the  84-day  treatment. 
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